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Abstract. Theoretical relative ages of galactic globular 
clusters have recently been challenged by a semi-empirical 
relation. It was used to point out that tested sets of iso- 
chrones were unable to reproduce the relation and yield 
internally inconsistent relative ages. We find that differen- 
tial cluster ages derived with the isochrones by Salaris & 
Weiss (1998) are reliable and internally consistent. We also 
show that this consistency depends on using the lower ab- 
solute ages determined by SW98, which therefore receive 
more empirical support. Moreover, we discuss the effect of 
the clusters absolute age on the evaluation of their differ- 
ential ages, and its connection with the question of their 
age dispersion. 
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1. Two opponent views of relative cluster ages 

The determination of galactic globular cluster (GC) ages 
has recently come into a state of flux, first due to up- 
dated stellar physics and lately due to the HIPPARCOS 
results. Theoretical isochrones were challenged by Buo- 
nanno et al. (1998, hereafter BCP), who have presented 
a self-consistent method for the determination of relative 
cluster ages, which as much as possible makes use of obser- 
vational properties and tries to minimize the input from 
theoretical isochrones. They argue that three tested sets 
of isochrones fail to reproduce the empirical relation be- 
tween A{B-V)^ B (the (B-V) difference between the TO 
and the base of the RGB) and [Fe/H] at a given age. Fur- 
thermore, relative ages based on brightness or colour dif- 
ferences are inconsistent in all cases. In this Letter we will 
confront the isochrones by Salaris & Weiss (1997, 1998; 
hereafter SW97 & SW98) with the observational results 
and show that this set passes the test to a great extent. 
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SW97 have determined GC ages by means of the tra- 
ditional theoretical approach of comparing theoretical iso- 
chrones with an observed CMD. The differences with re- 
spect to related works are (i) that they use the very lat- 
est (canonical) input physics, (ii) that also a-element en- 
hancement in the GC composition is taken into account 
(e.g. in the opacity tables) and (iii) that they determined 
absolute ages only for the very few clusters where both 
CMD morphology and available data guarantee an accu- 
rate dating. The brightness difference between Turn-Off 
(TO) and Zero- Age Horizontal Branch (ZAHB), AYgg , 
was used as the absolute age indicator. The imposed con- 
straints result in a rather small sample of 7 clusters of 
all metallicities (including disk GC in SW98) suitable for 
the determination of absolute ages (see SW97 and SW98 
for details). All other ages in their sample of 31 GC have 
been determined differentially with respect to these tem- 
plate clusters by using the (theoretical) dependence of 
A(B — V)^q B on age. For minimizing the errors aris- 
ing from uncertainties in the adopted colour transforma- 
tions and the mixing length calibration, SW97 and SW98 
evaluated the relative ages only within metallicity groups. 
Whether the relative ages would also be valid across metal- 
licity group boundaries, SW97 checked in only one case. 
This and other checks (e.g. availability of two clusters 
with an absolute age determination in the same metallic- 
ity group, comparison of the template clusters distances 
obtained from the theoretical ZAHB models with main- 
sequence fitting distances using HIPPARCOS subdwarfs) 
resulted in the confirmation of the internal consistency of 
the SW97 and SW98 ages and of the reliability of their 
theoretical isochrones. 

We now recall briefly the BCP method for determin- 
ing homogeneous relative ages for GC. The first step is 
to define a sample of coeval clusters in a wide range of 
metallicities. This requires age determinations, which at 
least must be able to define correctly what is "coeval" . To 
this scope, BCP used a variant of the vertical method used 
by SW97, employing a point on the main-sequence 0.05 
mag redder than the TO. The sensitivity of this new ab- 
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solute age indicator (AF - 05 ) and its reliability have been 
investigated in great detail in BCP. In total, a sample of 
11 clusters covering [Fe/H] = — 2.2 • • • — 0.5 was defined, 
whose absolute ages agreed within ±1 Gyr. While the ab- 
solute age of this group depended on the set of isochrones 
used (D'Antona et al. 1997; VandenBerg 1996, private 
communication to BCP; Straniero & Chieffi 1996, p.c. to 
BCP), membership within the group was basically invari- 
ant. From this, the empirical relation of A(B — V)^q B vs. 
[Fe/H] (at the common age of the group of clusters) is ob- 
tained (Fig. 11 of BCP), which in turn can be compared 
to the predicted one. Here, BCP noted that none of the 
adopted isochrone sets agreed with the empirical relation. 
It was also demonstrated that the discrepancy is inde- 
pendent of the colour-transformations employed and thus 
the inconsistency could arise from the theoretical models 
themselves. 

In the second step 6 clusters of definitely lower age 
than that of the calibrating group have been used to deter- 
mine the change in A(B — V) B q B with age at given [Fe/H]. 
The ages of the young clusters relative to the calibrat- 
ing ones were obtained by means of the Al/ 05 method. 
Again, theoretical models have to be used for this step; 
only theoretical predictions for the stellar brightness, but 
not for colours, are used. This finally yielded an empirical 
expression relating age differences with two observational 
quantities, namely A(B — V)^q B and metallicity (Eq. 3 
and Fig. 16 of BCP). This relationship is then applied for 
deriving relative ages for GC whose ZAHB level cannot 
be properly defined (blue or scarcely populated HB). For 
all other clusters in the BCP sample the average of the 
A t/0.05 and the A (£ _ t/)RGB relative ages is used. 

Considering that the BCP relationship depends itself 
on theoretical isochrones which were demonstrated to be 
inconsistent in at least one respect, it is worthwhile to 
question the reliability of the BCP-relation, which can be 
valid only if the theoretical isochrones do not suffer from 
systematic errors in the definition of the coeval clusters. 
Second, since SW97 and SW98 claim the internal con- 
sistency of their results, both in terms of absolute and 
relative ages, one could ask whether this result depends 
on the smaller sample of clusters selected by SW97 and 
SW98. In the following we will therefore re-derive the BCP 
relationship using the SW98 isochrones where theoretical 
input is necessary. 

2. Data set and coeval sample 

The reference GC data set is that of BCP with the ad- 
dition of NGC 6171 (a metal-rich halo GC included in 
SW97) and the exclusion of the three disk GC 47 Tuc, 
M 71 and NGC 6352. As extensively discussed in Alonso 
et al. (1997), A(B - V)%% B and AVgg yield significantly 
different ages for these clusters, if the typical GC he- 
lium content (0.24) is assumed. The problem can be re- 
solved with a nearly solar value (0.27; SW98) for the disk 
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Fig. 1. AV ' 05 for the 12 coeval clusters compared to the 
theoretical values by SW98 for several ages (in Gyr) 

clusters, but this would add an unwanted additional de- 
gree of freedom for the semi-empirical relationship investi- 
gated here. Moreover, given the strong dependence of the 
horizontal-branch luminosity on the He abundance, the 
inclusion of these clusters would affect the present com- 
parison with the models because it would affect the selec- 
tion of the sample of coeval clusters. The reader should be 
well aware of this difference between BCP and the present 
paper, because at the metal-rich end of the relationship 
A(B - t/)J r g B us[Fe/H] BCP found the largest discrepan- 
cies between observations and models. 

Adopting this data set we have selected a sample of 
coeval clusters, on the base of the observed AV om values 
and their errors, and the absolute ages as determined from 
the SW98 isochrones. The largest possible sample of coeval 
clusters is different from BCP and it is shown in Fig. [l] (see 
also data in Tab. |l]); individual ages range from 10 to 12 
Gyr and overlap within the errors, providing an average 
age of 10.9 Gyr. The sample spans the metallicity range 
-2.0< [Fe/H] <-1.0, typical of the bulk of halo GC, and, 
contrarily to BCP, it does not include any of the most 
metal poor clusters with [Fe/H]<-2.0. This is in agreement 
with the results by SW97 and SW98 who found that the 
most metal poor GC are on average older than the other 
GC. 

The result that the coeval clusters are not the same as 
in BCP is at odds with their conclusion that the member- 
ship to the coeval group does not depend on the isochrones 
considered. The main reason for this occurrence is the 
lower absolute age found for the individual GC with the 
SW98 isochrones. As previously described, the average 
age of the coeval GC is 10.9 Gyr, while BCP considered 
an average age of 15 Gyr. As it is evident from Figs. 7 
and 8 in BCP, the shape of the theoretical relation be- 
tween AT^ 05 and [Fe/H] depends on the absolute age of 
the isochrones. For example, at ages around 12 Gyr, the 
derivative 5AV 05 /5(a,ge) at low metallicities is quite dif- 
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Table 1. Data for the coeval GC sample 



Cluster 


AV ' 05 


A(B - V)?g B 


[Fe/H] 


NGC 362 


4 36+0 07 


263+0 01 


-1.27 


NGC 1261 


4.37±0.07 


0.261+0.01 


-1.29 


NGC 3201 


4.28±0.09 


0.243+0.01 


-1.56 


NGC 5272 (M 3) 


4.33±0.06 


0.249+0.01 


-1.66 


NGC 5904 (M 5) 


4.33±0.07 


0.250+0.01 


-1.40 


NGC 6101 


4.31±0.08 


0.239+0.01 


-1.81 


NGC 6121 (M 4) 


4.43±0.08 


0.248+0.01 


-1.36 


NGC 6171 


4.37±0.08 


0.260+0.01 


-1.00 


NGC 6362 


4.34±0.07 


0.264+0.01 


-1.08 


NGC 6584 


4.35±0.08 


0.248+0.01 


-1.54 


Pal 5 


4.27±0.08 


0.266+0.01 


-1.47 


Arp 2 


4.22±0.08 


0.248+0.01 


-1.84 



ferent from the case of 14 or 16 Gyr, independently of the 
set of colour transformations and bolometric corrections 
used. 

If we would have shifted by w +0.2 mag the AV - 05 
values of the isochrones used in BCP (see their Fig. 8) in 
such a way that the average age of the coeval clusters is 
around 11 Gyr, we would have excluded the most metal 
poor ones from the coeval sample, thus recovering the re- 
sult obtained with the SW98 isochrones. 

Another effect of the lower absolute age of the GC 
is that now we find included in the coeval sample two 
clusters, Pal 5 and Arp 2, considered as young by BCP. 

3. Differential ages 

The first step to evaluate GC relative ages following the 
technique by BCP is the calibration of a relation between 
A(B — V)to B an d metallicity for the subset of coeval 
clusters. A linear fit to the data (taking into account the 
A(B — V)to B error given in Tab.0 and a typical error of 
±0.15 dex in metallicity) provides 



A(B- 



t/\RGB 



(0.028±0.013)[Fe/H] + (0.294±0.020)(l) 



The next step is the evaluation of the ratio S/Atg, 
where S is the difference between the observed A(B — 
V)tc, B an( i the one expected on the basis of Eq. (1), and 
Atg is the relative age (in Gyr) with respect to the co- 
eval clusters, as determined from AV" ' 05 . We determined 
this ratio for the four younger clusters NGC 1851, Pal 12, 
Rup 106 and Ter 7, (see BCP for a discussion), then aver- 
aged these four values (having checked that they are not 
significantly correlated with metallicity) and get: 

S/At 9 = -0.0158 ± 0.0057 mag/Gyr (2) 

This quantity is smaller than that of BCP (-0.0093), im- 
plying a larger sensitivity of A(B — V)^o B to age and 
therefore smaller age differences for an observed A(B — 
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Fig. 2. Predictions of the SW98 isochrones (dashed lines ) 
for A(B — V)!^q B as a function of [Fe/H], compared with 
Eq. (1) (solid). Filled circles correspond to the observa- 
tional points for the coeval clusters. A typical error bar is 
shown as well) 



V)\^- By combining Eq. (1) and Eq. (2) we obtain the 
final relation between A[B - V)^% B , [Fe/H] and At g : 

At 9 = (—0.0158 ± 0.0057) -1 • [(A(B — V)^q B 

- (0.028 ± 0.013) [Fe/H] - (0.294 ± 0.020)] (3a) 

This equation allows the proper calculation of the error in 
Atg by standard error propagation. Taking into account 
the errors in A{B — V)t§ B an d [F e /H] as well, one gets 
errors of the order of ±2.0 Gyr in the derived values of 
Atg (see the vertical error bars in Fig. ||). Without the 
errors Eq. (3a) can be simplified to 



Atg = -63.3A(B - V) 



RGB 
TO 



1.8 [Fe/H] + 18.6 



(36) 



The coefficients depend of course on the absolute age of 
the coeval sample; for the sake of comparison we recall 
here that the corresponding values obtained by BCP were 
(-107.5, 4.3, 33.3). 

4. Discussion 

In Sect. 1 we recalled that BCP found an inconsistency 
in the theoretical isochrones by comparing the empirical 
relation Eq. (1) with them. In Fig. || we show the same 
comparison and, in the considered range of metallicity 
we did not find the same problem: the empirical relation 
crosses a 2-Gyr range in the isochrones, between 10 and 12 
Gyr, which is identical to the variation of ages determined 
from Al/ 0,05 (Fig. |l|). It is also evident that the spread of 
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-6.0 -4.0 -2.0 0.0 2.0 4.0 

At(SW98) 



Fig. 3. Comparison of the GC differential ages (in Gyr; 
with respect to t=10.9) as obtained by SW98 and from 
Eq. (3a). The dashed line displays the 1:1 relation, the 
solid one the regression line of the points. 



the observational points and the errors associated to them 
(see Table 1) is large enough to make the 11 Gyr isochrone 
fully compatible with the empirical relation. 

To further substantiate our claim that the SW98 iso- 
chrones - and therefore the semi-empirical relation de- 
rived in the previous section - are trustworthy, we show 
in Fig. H the comparison between the differential ages 
(with respect to t=10.9 Gyr) obtained with both methods 
for the 21 clusters in common (the coeval GC with the 
exception of NGC6121 and NGC6362, the 4 young ones 
and the blue HB clusters NGC288, NGC1904, NGC2298, 
NGC6254, NGC6397, NGC6752 and NGC7492) with 
-2.0 < [Fe/H] < -1.0 on the Zinn & West (1984) metal- 
licity scale, as adopted in the data set. Before doing so, we 
had to correct the ages found by SW98 for the different 
metallicity scale they have used, namely the Carretta & 
Gratton (1997) one. This new metallicity scale lowers the 
absolute ages by «0.8 Gyr when used instead of the Zinn 
& West (1984) one (SW98). Therefore, we have increased 
by this amount the ages given in SW98, and we have com- 
puted their differential ages with respect to the zero-point 
of the semi-empirical relative ages, fixed at 10.9 Gyr. 

The agreement between the two sets of differential ages 
is excellent: there is just a zero point shift of 0.4 ±0.5 Gyr, 
well within the errors associated with the age determina- 
tions. This small zero point shift reflects basically the fact 
that for the coeval clusters in common (10 out of 12), 
SW98 find an average age (after the correction for the dif- 
ferent metallicity scale) of 10.6 Gyr instead of 10.9 Gyr. 
When considering only the 7 blue HB clusters, one finds a 



zero point shift of 0.7 ± 0.8 Gyr, in good agreement with 
the results from the complete sample of 21 clusters. 

The dependence of Eqs. (3a-b) on the absolute age of 
the GC affects also the derived age spread among the clus- 
ters. To quantify whether a real intrinsic age dispersion 
exists among the clusters, over that due to the errors in 
the single age determinations, we have performed the F- 
test (Press et al. 1992) in the way described by Chaboyer 
et al. (1996) and SW97. It compares the observed age 
distribution with the expected one under the assumption 
of no intrinsic age dispersion, and provides an estimate 
of the probability that the two distributions are actually 
the same, taking into account the errors on the individual 
ages. If this test indicates that the sample is not coeval, the 
size of the true age dispersion (er rC ai) among the clusters 
can be estimated according to cr? cal = (<rl hs — <J^ p ), where 
Cobs and cr CX p are, respectively, the la dispersion in the 
observed age distribution and in the expected one. When 
considering the calibration of the semi-empirical method 
given by our Eqs. (3a-b), we find that the complete sample 
of GC (with the exclusion of the clusters with [Fe/H]<- 
2.0, for which we could not calibrate the method) shows a 
real la age spread by 0.9 Gyr. Performing the same test 
with the calibration of the method given in BCP for the 
same sample of GC, one obtains a la spread by w 2 Gyr. 

To conclude, we find that both the theoretical and the 
semi-empirical approach yield the same relative ages for 
the cluster sample of BCP when using the isochrones by 
SW98. In particular, the A(B - V)^ B values as a func- 
tion of age and metallicity predicted by these models are 
fully consistent with the observations. We also showed how 
the empirical method and results deduced from it (e.g. the 
amount of the age spread among GC) depend on the ab- 
solute age of the coeval sample and pointed out that the 
definition of this sample also depends on the GC absolute 
ages. The consistency of relative ages derived from the 
isochrones (Figs, [j] and |2|) and the agreement with those 
from the empirical relation is achieved only for the lower 
ages determined in SW97 and SW98. This result there- 
fore further substantiates the lower ages for GC as found 
recently in several papers. 
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